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Abstract. 

Reheating is a transition era after the end of inflation, during which the inflaton is 
converted into the particles that populate the Universe at later times. No direct cosmological 
observables are normally traceable to this period of reheating. Indirect bounds can however be 
derived. One possibility is to consider cosmological evolution for observable CMB scales from 
the time of Hubble crossing to the present time. Depending upon the model, the duration 
and final temperature after reheating, as well as its equation of state, may be directly linked 
to inflationary observables. For single-field inflationary models, if we approximate reheating 
by a constant equation of state, one can derive relations between the reheating duration (or 
final temperature), its equation of state parameter, and the scalar power spectrum amplitude 
and spectral index. While this is a simple approximation, by restricting the equation of state 
to lie within a broad physically allowed range, one can in turn bracket an allowed range of n s 
and r for these models. The added constraints can help break degeneracies between inflation 
models that otherwise overlap in their predictions for n s and r. 
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1 Introduction 

The inflationary paradigm [1-10] offers, in its numerous constructions (see e.g. [11]), a 
testable [12, 13] description for the physics of the very early Universe. Inflation addresses 
several open problems in cosmology, chief among them the question of the origin of cos¬ 
mological structures. In its simplest realization, the Universe is dominated by the potential 
energy of a light scalar field, the inflaton, that drives the expansion. In this picture, quantum 
fluctuations of the scalar field during inflation are precisely the primary source of cosmolog¬ 
ical perturbations [14-19]. The statistical properties of the Cosmic Microwave Background 
(CMB) fluctuations and of the Large Scale Structures (LSS) may therefore contain infor¬ 
mation about the physics of inflation. In addition to scalar density perturbations, inflation 
generically produces tensor perturbations, resulting in a spectrum of primordial gravitational 
waves which, via their impacts on the CMB and other astronomical sources, reveal informa¬ 
tion about inflation [20-23]. 

The transition from inflation to later stages of the evolution of the Universe (radia¬ 
tion and matter dominance) is referred to as reheating. During reheating the inflaton field 
loses its energy, eventually leading to the production of ordinary matter. Several reheating 
models have been proposed: the simplest ones, involve the perturbative decay of an oscil¬ 
lating inflaton field at the end of inflation [24-26], while more intricate scenarios include 
non-perturbative processes such as (broad) parametric resonance decay [27-29], tachyonic 
instability [30-35], and instant preheating [36 ] 1 . The word preheating indicates the initial 
stage of reheating, especially in the context where decay happens exponentially, generat¬ 
ing high occupation numbers in select frequency bands. Immediately after preheating the 
frequency bands that underwent parametric resonance will have extremely high occupation 

1 See also [37-40] for reviews and, e.g., [41, 42] for more studies on reheating. 


- 1 - 





numbers while the rest of the space will be basically un-populated, a highly non-thermal 
state. Over time, scattering events will spread out the distribution, eventually leading to a 
blackbody spectrum characterized by a final temperature T re , which normally corresponds 
to the temperature at the beginning of the radiation-dominated era. 

For some inflationary scenarios and for given interactions between the inflaton field and 
other matter fields, numerical studies were performed to derive an effective equation of state 
(eos). The eos is parametrized by a function w re {t ) for the Universe during the various stages 
of reheating. As inflation ends, the eos parameter is equal to —1/3. Assuming a massive 
inflaton, very quickly the eos climbs to 0, the eos of a massive harmonic oscillator oscillating 
between potential dominance (eos of —1) and kinetic dominance (eos of 1). During this initial 
phase of reheating, the frequency of oscillations, characterized by the inflaton mass rri, will be 
larger than the expansion rate. It is therefore correct to approximate the eos of the inflaton 
as a constant of 0. This is the equation of state of the Universe at the beginning of reheating 
when the Universe is still dominated by the inflaton field. As the inflaton decays and the 
decay products compose an increasing percentage of the energy density of the Universe, the 
eos will increase from 0 to 1/3 at the start of radiation dominance. In [43] it was shown 
that for a simple chaotic inflation model and for a quartic g 2 (j) 2 x 2 interaction (</> being the 
inflation and % its decay product), the equation of state right after inflation, characterized by 
w re = 0, sharply, within a couple efolds, changes to w re ~ 0.2 — 0.3 already during preheat¬ 
ing, long before the system reaches thermal equilibrium 2 . The duration of preheating can 
therefore generally be regarded as “instantaneous” in comparison with the remaining stages 
of reheating. In cases like the ones described in [43] (see also [27]), w re may therefore be 
rightfully treated as a constant throughout the entire reheating era. 

Aside from its thermalization temperature, T re , and effective equation of state, w re , 
reheating is also characterized by its duration, which one may quantify in terms of e-foldings 
N re = In (a re /a en d), occurring between the time inflation ends, t en d , and the beginning of the 
radiation-dominated era, t re . 

The reheating era is a difficult one to constrain observationally: except for some non- 
conventional scenarios (e.g. [45-56]). In the absence of topological defects like monopoles or 
strings, the fluctuations produced during reheating remain sub-horizon and cannot leave an 
observable imprint at the level of the CMB or LSS. A lower bound is placed on the reheating 
temperature by primordial nucleosynthesis (BBN) Tbbn ~ lO~ 2 GeV [57] 3 ; the scale of infla¬ 
tion is merely bounded from above (the CMB B-modes recently measured by BICEP2 [59, 60] 
do not yet, unfortunately, point to an inflationary signal) and can be as large as ~ 10 16 GeU, 
leaving for T re an allowed range of many orders of magnitude . Aside for the production of 
metric fluctuations in the aforementioned scenarios, a variety of signatures (or lack thereof) 
relative to the production of primordial black holes [61-63], magnetic field [64-66], unwanted 
relics [67, 68] and also to mechanisms such as baryo-and leptogenesis [69-73] (and more, see 
[39] for an overview and for a full list of related references), may be traced back to specific 
preheating/reheating models. 


2 A physical system reaching an effective (macroscopic) state characterized by nearly constant ratio of 
pressure over energy density while it is, microscopically, still out-of-equilibrium (“pre-thermalization”) had 
been previously investigated in Minkowski spacetime in [44]. 

3 Smaller values may be assigned to the lower bound of the reheating temperature in models such as [58]. 
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Another possibility for extracting information about reheating is to consider the expan¬ 
sion history of the Universe between the time the observable CMB scales crossed outside the 
Hubble radius during inflation and the time they later re-entered, in such a way as to define 
a relation between inflationary and reheating parameters [74] 


In 


k 

ao Hq 


—N k — N re 


Nrd + hi 


O’eqHeq 

do Hq 


+ In 



( 1 . 1 ) 


In this equation, k can be chosen as the pivot scale for a specific experiment, N k is the num¬ 
ber of e-foldings between the exit time of the modes at this pivot during inflation and the 
end of inflation, N re and Nrr> respectively indicated the e-folds between the end of inflation 
and the end of reheating and between the end of reheating and the end of the radiation- 
dominated era. From (1.1) one realizes that from the CMB constraints on the primordial 
power spectrum (which would correspond to a prediction for N k ), for a given inflationary 
model one would be able to infer the sum of Nrr and N re . To solve for N re and Nrr, 
individually one needs more information. For reheating models that can be parametrized by 
a constant effective pressure to energy ratio w re , one can relate the density at the end of 
inflation to the density at the end of reheating, and then assuming conservation of entropy 
after reheating, to the temperature today. This way one obtains another equation with the 
same two unknowns N re and Nrd that can be used to solve for each individually, or to 
rework the equations to trade the quantity N re for T re , the temperature at the end of re¬ 
heating. All of this is particularly straightforward for single-field models of inflation that are 
entirely defined by the form of their potential. In summary, for a given inflationary model 
and for given equations of state during reheating lying within a reasonable physically plau¬ 
sible range, one may use the CMB data to place constraints on the reheating temperature 
and its duration. These techniques have been successfully employed in several studies [75-82], 


In the same spirit as [75-82], and using similar techniques as in [81] (where the attention 
was directed specifically to inflation with power-law potentials, V(cj)) ~ </>“), we consider the 
constraints imposed by reheating on popular single held inhationary scenarios. We derive 
predictions for the length of the reheating era, and the temperature at the end of reheating 
for each model, assuming a constant equation of state during reheating. Accounting for the 
lower bounds on T re imposed by BBN and considering a physically plausible range of values 
for w re (likely the average value will fall between 0 and 5 ) we use the relations between 
reheating and inhationary parameters and the constraints on the primordial power spectrum 
amplitude and tilt from Planck [12, 13] to provide new constraints on the parameter space 
in given inhationary models. This is a useful and relatively new tool for constraining and 
differentiating between inflation models. Models might overlap in predictions for n s and r, 
but not for the same w re . As the constraints on n s gets tighter, this will translate into an 
increasingly narrow allowed range for w re for a given inflation model, and so this technique 
of constraining models with reheating will be increasingly efficient in ruling out some models 
in favor of others. 


This work is organized as follows: in Sec. 2 we detail the derivation of the reheating 
duration and of the temperature at the end of reheating as a function of the spectral index, for 
canonical single-held inhationary models and for reheating scenarios that can be described 
in terms of a constant effective equation of state; in Sec. 3 we review the analysis of [81] 
for a power law potential and we discuss the constraints from reheating on the inhationary 


- 3 - 











parameters; in Secs. 4 through 7 we compute the relations between inflationary and reheating 
parameters in the Starobinsky, Higgs, natural and hilltop inflation models and we discuss the 
bounds placed on some of these models by reheating; in Sec. 8 we present our conclusions. 

2 Calculating N re and T re 

A reheating model (or class of models) may be characterized by a thermalization temperature 
T re , a duration, N re (here defined in terms of the number of e-folds counted from the end 
of inflation), and an equation of state with an effective pressure-to-energy-density ratio, w re . 
The latter should have values larger than —1/3 for inflation to come to an end, and is assumed 
to be smaller than 1 in order not to violate causality. A variety of reheating scenarios allow 
for an equation of state that is nearly constant in time. For the purposes of this work we will 
thus approximate w re as a constant in all our calculations; in our plots for N re and T re , we 
assign to w re sample values ranging in the interval [—1/3,1]. We define N re as the time frame 
from the end of inflation until the equation of state makes a step function transition from 
the value w re it had during reheating to w = 1/3, which we define as the start of radiation 
dominance. T re is the temperature when this transition occurs. From this definition, N re 
and T re are not well defined if the equation of state during reheating is also equal to 1/3 
(we will discuss this case more later). Also, we assume a standard expansion history after 
reheating, with a radiation-dominated (RD) era followed by a matter-dominated (MD) one. 

We derive, following [75-82], an expression for the reheating parameters ( N re , T re and 
w re ) in terms of a set of physical quantities that are specific to inflation and to the cosmolog¬ 
ical epochs subsequent to reheating. Considering the evolution of the Universe between the 
Hubble-exit time during inflation (henceforth indicated by tk) for observable scales and the 
time of observation of the same scales (to), one can write matching conditions for the total 
energy density as well as for the scale factor, a(t), during the intermediate eras. Fig. (1) 
summarizes the evolution of the comoving horizon distance throughout this length of time, 
marked by the transitions between consecutive epochs at t en d, the end of inflation, t re , the 
end of reheating/beginning of RD era, and t eq , the beginning of the MD era. 

In the figure we equate the size the comoving horizon far back into inflation, correspond¬ 
ing to modes l = 2, to the size of the horizon today. In order to solve the horizon problem, 
the span of comoving scales that leave the horizon from l = 2 to the end of inflation must 
equal the span of comoving scales that reenter the horizon after inflation till today. Note the 
factor by which the comoving horizon shrinks between scales l = 2 and the end of inflation 
(the length of the first line in the figures) is not known. The slope of that line is set by the 
fact that the equation of state is ~ — 1 during inflation. Depending on the model, that line 
could be longer or shorter. While there is a minimum length in order to solve the horizon 
problem while having Inflation occur before BBN, there is no upper bound. The value of 
w re will set the slope of the second line, the rate by which modes reenter the horizon during 
reheating. In the figure we display the two extreme cases of w re = 1 and w re = —1/3. One 
can see from comparing the two plots, the smaller w re is during reheating, the less efficiently 
modes re-enter the horizon, and the more efolds will be necessary in the post-inflation period. 

We consider single-field inflationary models with background field equations, cj) + 3H<j) + 
V = 0 and 3 H 2 Mp ~ V(4>). We also assume that both e and remain smaller than 1 
throughout the inflationary regime. 
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Figure 1: Each figure shows the evolution of the comoving horizon distance over time. 
Each figure shows the extreme cases for w re : the first figure for w re = 1 and the second for 

W r e — 3 ■ 


If one assumes a constant equation of state, the change in the scale factor during reheating is 
easily related to the change in the energy density. Using p oc the reheating epoch 

is described by 

Pend / ^end 

Pre \ C^re 

where the subscript end refers to the end of inflation (the start of reheating), and re refers 
to the end of reheating. Writing this in terms of e-foldings 


—3(l+«V e ) 


( 2 . 1 ) 


N re 


_1 - ln f Pend\ 

3(1 + Wre) \ Pre J 


3(1 -{- w re J \2 p re J 


( 2 . 2 ) 


where the last step of (2.2) is obtained by replacing p en d = (3/2 )V en d, derived by setting 
w = —1/3 at the end of inflation. 

The temperature is related to the density by 

7T 2 

Pre = 7^9reT re , (2-3) 

where g re is the number of relativistic species at the end of reheating. Combining Eqs. (2.2) 
and (2.3) one finds 


1 / 30 • | V end \ 

3(1 + w) n ( Tr 2 g re T? e J ' 


(2.4) 


Making the standard assumption that entropy is conserved between the end of reheating 
and today, one can relate the reheating temperature to the temperature today by taking 
into account the changing number of helicity states in the radiation gas as a function of 
temperature, 


T re — Tq 




(2.5) 
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where Nrd is the length in e-folds of radiation dominance, e Nrd = a re /a eq . The ratio 
a 0 /a eq can be rewritten as 


ao_ = apH k c - Nhc - N „ c - Ntin . 

CLeq k 


( 2 . 6 ) 


where one uses the relation k = a k H k for the time at which the pivot scale k 1 crosses outside 
the Hubble radius and N k is defined as the number of e-foldings between the latter and the 
time inflation ends. Inserting (2.6) into (2.5) one finds 


T 

- 1 re 



a 0 T 0 

~k~ 


H k e 


- N k e -N r . 


(2.7) 


Notice that larger values of N re corresponds to smaller T re and vice versa. In other words, as 
expected, the quicker and more efficiently reheating takes place, the larger the temperature. 
Plugging (2.7) into Eq. (2.4) 


N r - = 


3(1 + w re ) 


s ln 


3 2 • 5 

^~9re 


+ In 


YZd 

He 


1 . 

H— In 
3 


1 l*?n 
”13” 


+ In 


«o Tp 


+ N k + N r , 


( 2 . 8 ) 


One can first solve for N re assuming w re i 


N re = 


(1 3 W r e) 


"I 1 " 


3- -5 

^9re 


- - In 
3 V 43 


11 9r 


— In 


apTp 


— In 


H k 


~N k 


(2.9) 


Notice that the values of the last two terms in Eq. (2.9) depend on the specific inflationary 
model. Assuming g re ~ 100 and using Planck’s pivot of 0.0 bMpc^ 1 5 , one obtains a simplified 
expression for N re , before specifying a particular inflationary model: 


Nre = 


(1 3w r e ) 


61.6 - In | 

H k 


~N k 


( 2 . 10 ) 


One can then use Eq. (2.7) to obtain 


T 

- 1 - re 



5 


3 2 • 5V en d 

l i 

3(1 -\-w-re) 

TT^gre 



3 (1 —| -Wre) 
3w re — 1 


( 2 . 11 ) 


4 Note in the following when we repeatedly refer to the pivot scale, we will use throughout Planck’s pivot 
scale of 0.05Mpc _1 . 

s The convention in the Planck analysis defines the pivot scale such that the comoving momentum k becomes 
horizon sized when kag = aH, where we have been using k = aH, so using our conventions N = 0.05Mpc -1 . 
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2.1 Special case w re = | 

The final result for N re in Eq. (2.10) only applies for w re ^ 1/3. Going back to Eq. (2.8), 
notice that if w re = |, N re cancels from both sides of the equation, and one is left with 



K Z g r e J 


+ 7 In 
4 




H— In 
3 


f 11 9re 

V 43 



+ Nk . 


( 2 . 12 ) 


Assuming g re = 100, and Planck’s pivot scale, this simplifies to: 


61.6 = In 


H k 


+ AIfc . 


(2.13) 


For w = 1/3, it is not possible to derive a prediction for N re or T re but instead, for a 
particular inflation model, one finds a prediction for n s . Note the ambiguity in N re and T re 
is due to the fact that we are defining the start of radiation dominance as the moment w re 
reaches 1/3. If w re is already equal to 1/3 during reheating, then there is ambiguity in when 
to differentiate between the two regimes. 


2.2 Model dependent part 

In order to solve for N re in Eq. (2.10) (or to solve for n s in Eq. (2.13) if w re = 1/3) for a 
particular model, one needs to compute N k , H k , and V en( i- N k can be calculated starting 
from the definition of e-foldings: 


AN = 



(2.14) 


Recasting the r.h.s. of (2.14) as an integral over <fi and using the background equation of 
motion for the inflaton, 3 H<f> + V' ~ 0, and the Friedmann equation, H 2 ~ V/(3Mp), one 
finds 

1 y 

(215) 

Next, H k can be written as a function of n s . Using the definition of the tensor-to-scalar ratio 
r = Ph/P( (where P k = (2H 2 ) / (it 2 Mp) and P ,/ = A s at the pivot scale) 


r k 


ZHj 

TT 2 MpA s 


Then using r = 16e this gives 

H k ~ irM P ^8A s e k . 


(2.16) 


(2.17) 


Once the form of V(4>) is specified for a given model, one can express V en d as a function of 
model parameters calculated at the pivot scale. The explicit form of V en d along with (2.15) 
and (2.17) can be plugged into Eqs. (2.10) and (2.11) to derive N re and T re as a function of 
inflationary model parameters (or into Eq. (2.13) in the case w re = 1/3). 


-7 - 









3 Polynomial potentials 

Consider a polynomial type potential 


V = (3.1) 

This was considered in the context of reheating in [11, 75, 77, 81, 82]. We quickly review this 
specific application. At the end of this section, we discuss with some quantitative examples 
how closely the constraints from inflation compare to the ones from reheating. 

The first step is to calculate the model dependent parameters in Eq. (2.10), i.e. Nk, Hk, and 
V end - The number of e-folds between the time the pivot scale exited the Hubble radius and 
the end of inflation can be derived using Eq. (2.15) 

W ■ (3 - 2) 


The potential in these polynomial models is generally steep enough so that </>& (j> en d and 
it is appropriate to approximate 


N k 


2aMp 


4>l- 


(3.3) 


We now require Nk as a function of n s . From the expression of the spectral index as a function 
of the slow-roll parameters, n s = 1 — 6 e + 2 rj (where e = (Mj,/ 2)(V'/V) 2 and 77 = MpV"/V), 
and using (3.3) to rewrite e and 77 as functions of Nk, one finds 


N k 


a + 2 
2(1 -n s y 


(3.4) 


From Eq. (2.17) and using the previous equation, Hk is given by 


I A'jr 1 

H k =irMpJ—^(l-n s ). (3.5) 

Lastly one computes V en d in terms of n s and A s , 

V end = 3M 2 H 2 ^ = 6 tt 2 M p A s (1 - n s ) ^ ) , (3-6) 

where the value of the inflaton field at the end of inflation was computed by solving for 4> end 
from the condition e = 1 . 


Thus Nk, Hk, and V end are all expressed as functions only of a, n s and A s and one 
may plot N re (and T re ) as a function of n s for some fixed values of w re and a. We use 
n s = 0.9682 ± 0.0062 and Planck’s central value A s = 2.196 X 10 ~ 9 (small variations in A s 
have negligible effects on reheating predictions). 

We plot in Fig. 2 N re and T re predictions for a = 2/3, 1, 2 and 4. The case a = 2/3 is 
favored by axion-monodromy models, and a = 1 and a = 2 give promising predictions when 
compared with the Planck data. The case a = 4 is difficult to reconcile with w re < 1 even 
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Figure 2: Plots of N re and T re , the length of reheating and the temperature at the end 
of reheating respectively, for polynomial potentials with exponent a. The solid red line 
corresponds to w re = —1/3, the dashed green line to w re = 0, the dotted blue line to 
w re = 2/3, and the dot-dashed black line to w re = 1. The pink shaded region corresponds to 
the la bounds on n s from Planck. The purple shaded region corresponds to the la bounds of 
a further CMB experiment with sensitivity ±10 -3 [83, 84], using the same central n s value as 
Planck. Temperatures below the dark green shaded region are ruled out by BBN. The light 
green shaded region is below the electroweak scale, assumed 100 GeV for reference. This 
region is not disallowed but would be interesting in the context of baryogenesis. 

considering the 2a bounds on n s 6 . 

Instantaneous reheating is defined as the limit N re —> 0, visualized in the figure as the point 
where all the lines converge. Such instantaneous reheating leads to the maximum temperature 
at the end of reheating, and the equation of state parameter is irrelevant. 

(Thus, while not shown, a w re = ^ solution would correspond to a vertical line passing 
through the instantaneous reheat point.) 

From Fig. 2, a = 2/3 can be consistent with Planck bounds, but assuming an equation 
of state w re > 0, the model would tend to predict smaller reheating temperatures if one 
considers Planck’s la bound on n s \ using Planck’s 2cr bounds, any reheating temperature up 

6 An exception where </> 4 may still be viable is in the context of warm inflation [85, 86]. 
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to the maximum instantaneous case is still allowed. 

For a = 1 and a = 2 all the lines in Fig. 2 are shifted towards the central value of n s 
when compared to the a = 2/3 case, thus allowing for a wider range of reheating tempera¬ 
tures as well as values of the equation of state parameter. 


Consider now the case w re = 1/3. Solving Eq. (2.13) for the polynomial potential, one 
obtains 


61.6 


1 / 3a \ a + 2 

4 \^'K 2 A s {a + 2)) + 2(1 - n s ) 


(3.7) 


Using Planck’s central value for A s , Eq. (3.7) gives specific predictions for n s 


n s = 0.977 
n s = 0.974 
n, = 0.965 


for a = 3 , 
for a = 1 , 
for a = 2 . 


(3.8) 


Notice that larger values of a require smaller values of n s . With the 2a bounds on n s from 
Planck, 0.956 < n s < 0.981, w re = | would be consistent with all three values of a. 



Figure 3: Parameter space for 0 2//3 inflation. The figures show r and Nk predictions that 
give the correct A s for the plotted n s at the pivot scale. The green portion of the line 
comprises the region of parameter space corresponding to reheating models with w re > 1 , 
the yellow part corresponds to w re >1/3, red to w re <1/3 and orange to w re < 0. Note the 
most likely w re , between 0 and 1/3, falls in the red region. 



Figure 4: Parameter space for <f> inflation. Shading is as for Fig. (3). 















0.25 



Figure 5: Parameter space for 0 2 inflation. Shading is as for Fig. (3). 


Figs. 3-5 shows the parameter space in the r and N[. vs. n s plane, corresponding to 
the different reheating scenarios. We allow for any Nk > 19. We note again that there is 
no maximum allowed N^. A minimum on A/ is determined by the temperature at the end 
of reheating in order to solve the horizon and flatness problems. One finds that N > 24.9 if 
reheating after inflation is to be above the BBN scale and N > 34.8 for reheating above the 
electroweak scale in order for scales on the order of the horizon today ( i .el = 2) to have left 
the horizon during inflation. A simple estimate of the ratio of expansion scales between 1 = 2, 
and Planck’s pivot scale, at l ~ 685.8, if the expansion rate during Inflation were constant, 
is AN ~ Info/h) ~ 6.5. However, in the large field modes we are considering, the variation 
in H is not negligible and the exact AN = ln( ^j^ ) is closer to AN ~ 5.9. This means that 
for reheating greater than the BBN scale, one finds Nk > 19 (or Nf. > 29 for reheating above 
the electroweak symmetry breaking scale). 

The green part of the line in Fig. 3 corresponds to the region of parameter space that 
requires reheating models with w re larger than one, the yellow part corresponds to w re >1/3, 
red to w re <1/3 and orange to w re < 0. We stress that a value of w re between 0 and 1/3 is 
most likely and these solutions fall in the red band in Fig. 3. 

One can see that requiring 0 < w re <1/3 corresponds to respectively setting an upper and 
a lower bound on the tensor-to-scalar ratio 

0.05 <r< 0.06 for a = §, 

< 0.07 < r < 0.09 for a = 1, (3.9) 

0.14 < r < 0.18 for a = 2. 

Since it now appears that the majority of BICEP2’s signal is comprised of dust [87], it is 
difficult to find a viable reheating scenario for 0 2 inflation; if we loosen our restriction to just 
requiring w re < 1 then one obtains a bound r > 0.11, which is just inside the 2 a limit [87]. 

The assumption 0 < w re < 1/3 results in tighter constraints on r than Planck’s 2a 
bound on n s alone. For cf> 2 , the n s 2a bound yields 0.08 < r < 0.18. Restricting w re also 
provides stronger constraints on A>: for q i 2 , the n s 2a bound yields 45 < Nk < 103, whereas 
0 < w re <1/3 yields 44 < Nk < 57. 
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4 Starobinsky model 


The action for the Starobinsky model [1] has the form 


S 



'Ml 

2 


(R T OiR ) + £matter i 


(4.1) 


where R is the Ricci scalar. Performing a conformal transformation [88, 89] 

9fn/ = ^ (4.2) 

where cu 2 = 1 + 2 aR, the action (4.1) is rewritten as the canonical Einstein-Hilbert action 
plus other terms which form a modified C ma tter 


S 



~ M P [ s _ m ft 

2 [ (1 + 2 act)) 2 


6a 2 

(1 + 2 a^) 2 



+ £ 


matter ? 


(4.3) 


where what we now call <f> is equal to R, the original, untransformed Ricci scalar. Notice that 
d a carries factors of the metric, therefore ^ d a . Next one defines cj), a canonically normalized 
version of 


(j> = 



Mp ln(l + 2 acj)). 


(4.4) 


Rewriting the action in terms of </> one finds 


S = 


d A x 


’ Ml 


1 ( 

[2 0 \ 2 

1- x Q _ 0, /2~ 4> 

~Y 

R- 

Ta V 1 " 6 

/ 3 M p ] 

- 2 (^^) 2 + e ^ 3 Mp Lmatter 


• (4-5) 


If one assumes that the other fields in £ ma tter are subdominant during inflation and can be 
ignored, then one can verify that this Einstein frame action behaves as normal gravity plus 
a canonical scalar field with the potential 



(4.6) 


Dropping the bar on cj) from now on, but continuing to work with the canonical version of 
the field, one can easily compute the number of e-foldings between the horizon exit of the 
pivot scale and the end of inflation 


N k 


1 

Ml 





4>k 

(fiend, 


(4.7) 


/2 th 

With the approximations 4>k 3 > 4> e nd > an( 4 Mpe V 3 m p cf> k , the previous expression simpli¬ 
fies to 


which can be inverted for (j) k 



4>k 



Mp In 



(4.8) 


(4.9) 
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The next step is to compute e*, and rj k in order to derive N k as a function of n s using 
n s = 1 — 6e + 2r/. The slow-roll parameters have the following form 


tfc - 


4iV|’ 


Vk = 


X’ 


(4.10) 


where Eq. (4.9) was used along with the approximation N k 1. From Eq. (4.10) then one 
finds 


N k = 


1 ~ n 8 

Using the expressions above, one derives H k as a function of n s and A s 


H k = irMp\l -A s (l - n s ), 


V end = ^Tr 2 Mf,A s (l-ns)' 2 

2 (1-Ki-n,)) 2 


(4.11) 


(4.12) 


(4.13) 


Eqs. (4.11)-(4.13) are all that is needed to derive the results for the duration and for the 
temperature of reheating. 



Figure 6: Plots of N re and T re , the length of reheating and the temperature at the end of 
reheating respectively, for Starobinsky and Higgs inflation. All curves and shaded regions 
are as for Fig. 2 

Fig. 6 shows good compatibility with Planck’s lcr bounds on n s for all the possible w re 
values. Also, if one does not put any restrictions on the value of w re then any temperature 
between the BBN bound and the instantaneous reheating value is allowed within the la 
bound. 
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Figure 7: Parameter space for Starobinsky inflation. Shading is as for Fig. 3. 

We plot in Fig. 7 the parameter space in the r and TV/. vs. n s plane for Starobinsky inflation, 
for different ranges of w re ". For 0 < w re < 1/3, the corresponding range for the spectral 
index is 0.953 < n s < 0.964. This also corresponds to the range 0.004 < r < 0.007 and 
42 < N k < 56. 


5 Higgs Inflation 


The idea behind Higgs inflation [90] is to allow the Standard Model Higgs held to be the 
inhaton by adding a non-minimal coupling to gravity. The Jordan frame action is 


S 




1 + 2£ 


H\ 

~Ml) 


+ C 


matter > 


(5.1) 


where H is the Higgs doublet. We may again perform a conformal transformation to write the 
action in the form of Einstein gravity plus a modified Cmatter■ The transformation is given 
by g^u = u 2 g IJM with w 2 = 1 + 2Rewriting the action in terms of the transformed 
metric, we find 


5 



[ M Pp 3£ 2 

_ 2 u 4 Mp 


[dH^Hy 


+ - a C'matter 
Cc/ 4 


(5.2) 


Next, one extracts the kinetic and potential terms for the Higgs held contained within C ma tter- 
One can use V/ = H — i y) 2 , dropping the v part (we are interested in inflation scales 
much larger than electroweak scale). Ignoring all the Higgs interactions with other helds, 
and only considering its self coupling (which we assume is the dominant term in the Higgs 
potential at inflation scales) 


S = J d 4 xy/-^ 


—- 


3£ 2 


u 4 Mp 




Wh) +—C, 




matter 


(5.3) 


where now Cmatter comprises all the matter helds except the Higgs. Note: one needs to 
convert d a —> ui 2 d a . The Higgs is no longer canonical because of the effect of the non- 
minimal coupling. To canonically normalize all four of the Higgs degrees of freedom, one 

7 See also [11, 82] for n s vs r plots in the Starobinsky model for w re = 0. 
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must work in unitary gauge, where three of the four degrees of freedom are equal to zero, 


H = A= 


0 


V2 V h , 

S = I d 4 x\/—i 


M 2 p ~ 3 i 2 h 2 

2 K w 4 M 2 


dh ) 2 - -^(dhY - -A/i 4 + ^-C 


2oj 2 


4 uj 4 




matter 


(5.4) 


The canonically normalized version of h is h, defined as 


dh 

dh 




M°p 


h 2 (6£ + 1). 


(5.5) 


Before integrating the previous equation, it is useful to introduce a few approximations. First 
one uses 6£ 1. To get a successful inflation model, one should require £ ~ 10 4 . Next one 

uses the condition (6£)/(Mj,)/i 2 1. h ~ Alp when inflation ends, and therefore h > Alp 

for the duration of inflation. This allows one to rewrite Eq. (5.5) as 


h = ^f dh - h 


M P 


l + £il 
1 + Ml 


which integrates to 


(5.6) 


h= \j l M P In ( 1 + ^2 


(5.7) 


Rewriting the action in terms of h, one finds 


S= d 4 xy/-I 


^ A - l -(Sh\ 2 - XMi 


4« 2 


1 — e 


— 2 


P ) +e V 3 M p C-matter 


, (5.8) 


The potential term for the canonical held takes the same form as the Starobinsky potential 
with the identification (l)/(8a) = (AMj,)(4£ 2 ). Since we have a canonical held evolving 
in the same potential as the Starobinsky case, the Higgs inflation model gives the same 
predictions for N re and T re (see also [11, 82]). We note Starobinsky and Higgs inhation have 
different low scale behavior [91-95] and so while the allowed parameter space as a function 
of w re is the same, the w re that is most likely for Starobinsky vs. Higgs inhation is likely to 
differ. Tighter constraints could be obtained by considering gravitational, Planck suppressed 
couplings in the Starobinsky case, and standard model couplings in the Higgs case [91-93]. 
Of course new physics may modify the running of the couplings or add new couplings at 
these high scales (see for example [94]); in this respect, our approach of characterizing an 
allowed parameter space by assuming a range of w re between 0 and 1/3 can usefully help 
bracket different allowed scenarios. 


6 Natural Inflation 


The potential for natural inhation is [96] 


V{<f>) = A 4 


1 + cos 



( 6 . 1 ) 


- 15 - 














The number of e-folds N k between the time the pivot scale modes crossed outside the horizon 
and the end of inflation is given by 


N k = 


{ J_ V ln [~ sin 2 {Xend/2) 


\M P J L s in 2 (Xm/2) 
where x = </>//• The slow-roll parameters have the following form 


e = 


1 /M P N 


1 - cos(x)" 

2 V / z 

I 

1 + cos(x)_ 


V = 


f Mp\ 2 cosx 
\ f J 1 + cos x' 


( 6 . 2 ) 


(6.3) 


The field value at the end of inflation can be determined by setting e = 1; this leads to the 
following equation for Xend 


1 Mj, sin 2 ( Xend ) 

2 f 2 [1 + COS ( Xend)Y 


= 1. 


The solution is 


COS (Xend) — 


-1 + 6 

1 + 6 


\V2f) ' 


The number of e-folds in Eq. (6.2) can be written as 


N k = 


_ ( / Vi r 1 - cos(xend) ] _ ( f \ 2 


In 


1 


\MpJ L 1-COS (Xm)J \MpJ [(i + b ) i 1 - cos (Xrn)) 

The value of the field at the pivot scale during inflation is then given by 


z = 


COs(\'i n ) = 1 - Z , 

Using (6.3) and (6.7), one finds 

n s — 1 = —6e + 2r] = 


(1 + 6 ) 


exp 




(6.4) 


(6.5) 


( 6 . 6 ) 


(6.7) 


_ / Mp \ 2 (2 + 


V / J V 2 -zj ’ 


which leads to 


N k = 


-(A) 


2 

1 In 


1 + 


Mr 


(1 


Mp 

f 2 . 


2 /v l(i - ns ) + 


( 6 . 8 ) 


(6.9) 


Notice that the previous expression is positive and real only if the argument of the logarithm 
is defined between zero and one 


0 < 1 + 


Mf\ / (l ~n a )~7 


Mj, 


V 2 ) 


(1 — n s ) + -j 


Mj, 


< 1. 


( 6 . 10 ) 


The conditions (6.10) are equivalent to requiring that 




> 


i 


\M P J (1 - n s ) 


and 


3 + n s + 




( 6 . 11 ) 


16 



























The second condition in (6.11) is always true. The first condition implies a minimum / for 
each n s , and the bound on / increases with increasing n s . Using the central value for the 
Planck constraints on the spectral index, then (6.11) gives / > 5.6 Mp. 


/ ~i n. Mp 

(! -n,) - -y- 


The tensor-to-scalar ratio can be expressed in terms of n 

r = 16e = 4 

Using Vin = A 4 [1 + cos(xm)] — 3H' 2 Mp, one finds 

i + (i-0 (sfc) 


( 6 . 12 ) 


V end = 3H 2 Mp 
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2 + 4 (m?) 


(6.13) 
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Figure 8: We show N re and T re , the length of reheating and the temperature at the end of 
reheating respectively, for natural inflation, for 3 values of the coupling /. Again, shading is 
as in Fig. 2 

Fig. 8 shows N re and T re solutions for various reheating parameters w re and for various 
couplings / in natural inflation. Unlike polynomial inflation, or Starobinsky/ Higgs inflation, 
natural inflation has an extra free parameter, and so one no longer gets a precise prediction 
for the temperature and length of reheating once a reheating model, w re , and n s are specified. 
But one can get reasonable bounds on the coupling / such that a viable reheating model 
exists. 

One can obtain separate, and stronger constraints on / based on the requirement for 
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viable reheating. These constraints likewise are functions of n s , and their effects are displayed 
in Figure 9 8 . 

There is no upper limit on /. For / > 14 Mp the various w re lines reach an asymptotic 
form. As a result even for very large /, there is a valid solution for each w re value consistent 
with Planck’s la bounds. The asymptotic solution for large / for 0 < w re <1/3 corresponds 
to a solution for the spectral index in the range 0.956 < n s < 0.965. 



0.00 1 — ‘ — i —*—■■ lMl .. i ■ i i i i i m i - i ol ‘ - ■ ■ ■ 1 ■ ■ ■ ■ ‘ ■ ■ ■ ■ ‘ ■ ■ ■ ■ ‘ ■ ■ ■ ■ i ■ ■ ■ ■ i ■ ■ ■ ■ 1 

0.93 0.94 0.95 0.96 0.97 0.98 0.99 1.00 0.93 0.94 0.95 0.96 0.97 0.98 0.99 1.00 

n s n s 

Figure 9: Representation of the r and Nk vs. n s plane for natural inflation. The colored 
region is the entire allowed parameter space that can produce the measured A s value at the 
pivot scale, and the n s value at the pivot scale as plotted. Shading again follows Fig. 3 

Note the minimum on / increases with increasing n s , such that for larger n s , a larger / 
is needed to find a solution consistent with the reheating model being considered. The top 
limit in Fig. 9 A (the bottom limit in Fig. 9 B) is approached asymptotically for large /. The 
bottom part of the parameter space in Fig. 9 A (top part in Fig. 9 B) corresponds to small 
/. Everywhere in the figure > 19, such that inflation lasts long enough to allow for BBN. 

Using Planck’s 2a bounds on n s , requiring w re < 1 gives r > 0.02, and requiring 
w re <1/3 gives r > 0.05. For w re < 1/3, values of n s smaller than Planck’s central value 
would be favored for any /. The weakest constraint is for large values /, for which n s < 0.965. 


7 Hilltop inflation 

The potential is given by [4, 97] 


V (<(>) = M 4 


1- I - 

h 


pi 


(7.1) 


We begin by considering p > 2. The exact expression for the number of e-foldings 
between the time the pivot scale crossed outside the horizon and the end of inflation is 


N k = 




2 pMl 


P L 


2 2 i 

Xin Xend < 


P~ 2 


2—p 

Xin 


p — 2 


2—p 

Xend 


(7.2) 


3 See also [11, 82] for n 3 vs r plots in the natural inflation model for w re = 0. 
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where one defines x = ^/p- 
The slow-roll parameters are 


Setting e 


P 2 M'p X 


e = 


2(p—1) 


_ / i Mp 

V = ~p(:P ~ 1 )' 


X 


p -2 


2 p 2 ( l-x p f 1 V 2 (! — X p )' 

1 at the end of inflation one derives the equation for Xend 


p 2 A/J 


2(p-l) 

Xend 


^ (i-XendV 


= 1. 


(7.3) 


(7.4) 


Let us consider the case where p > Mp and define q = Mp/p. For small values of q, one can 
search for a solution for Xend hr the form of a Taylor expansion around q = 0 


Xend ®0 T 0-1 q T 


l -a 2 q 2 + 0{q*). 


(7.5) 


One can show that, up to order q 2 , a solution to Eq. (7.4) is (see e.g. [11]) 


1 „ , iP~ !) „2 

Xend — 1 9 T ^ q ■ 


(7.6) 


Similarly, one can look for a solution for the initial value of the scalar field using (7.6) and 
(7.2), to find 


Xin — 1 


l + 4JV fc 


q + 0{q 2 


(7.7) 


Using (7.6) in the expression for the slow-roll parameters, (7.3), the spectral index as a 
function of is 


n, — 1 ~ —; 


6 


1 + 4AV 


N k ~ 


1 

4 


f— 

V 1 - n. 



(7.8) 


The tensor-to-scalar ratio is 


8 n r 
r ~ -(1 - n s ). 

Notice that (7.8) and(7.9) only apply for small values of q, more precisely for 

y/2 

Q (p-l)Vl + 4iV fc ' 

For p G (3,8) and for N k G (30,100), the previous condition is satisfied if q < 0.01. 
Within the same range of validity, the potential at the end of inflation is given by 


(7.9) 


(7.10) 


Vend — 


\j\- h 2 m 2 VT= 


n,. 


(7.11) 


For p < 2 one derives the same results as Eqs. (7.8)-(7.9) and (7.11). For p = 2, however, a 
new expression for N k is required. In this case we find 


N k 


h 2 \xL 

2 M 2 [ 2 


Xend 


2 


In Xin + In Xend 


(7.12) 
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Figure 10: Plots of N re and T re , for hilltop inflation with p = 2 and for three different 
values of p. Shading is as in Fig. 2 

The plots for hilltop inflation are derived using a numerical procedure and therefore they 
convey more information than one would obtain with the above analytic results since they 
cover a range in which the latter would not apply (i.e. for smaller values of p) 9 . 

Note that for p = 1 the potential is just a straight line, and so should give the same predic¬ 
tions as the V (x (j) inflation model considered above. In Figs. 10-12 we therefore plot N re 
and T re for various reheating scenarios parametrized by w re , and various values of p for p = 2 
and larger. 

Just as with natural inflation, hilltop inflation has two free parameters, in this case M and 
p. This extra freedom means for each different p value, there are p values that are readily 
consistent with Planck data and p values that are not. One can give bounds on p for each p 
model such that reasonable reheating solutions exist, and these results are shown in Figure 
13 and 14. 

As with the bound on / for natural inflation, there is a minimum on p required for p = 2 
to get any solution at all, even before reheating is considered, and that bound is a function 

9 Notice that hilltop inflation was previously studied in the context of reheating in [11, 77, 82]. 
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Figure 11: Plots of N re and T re , for hilltop inflation with p = 3 and for three different 
values of p. Shading is as in Fig. 2 


of n s . This p m i n corresponds to r —> 0, so in this case there is no minimum on r before 
reheating is considered. For p = 3,4 there is no such minimum on p to get a solution in the 
regime p> Mp. There appears to be no observational constraint from Planck for very large 
values of p. 

The upper bounds in the plots in Figure 13 (the lower bounds in Figure 14) correspond 
to larger p, and the lower bounds in Figure 13 (or the upper bounds in Figure 14) correspond 
to smaller p. Using the 2<r bounds on n s and requiring w re <1/3 gives the following lower 
bounds on the tensor-to-scalar ratio: r > 0.02 (for p = 2) and r > 0.007 (p = 3) and r > 0.003 
(p = 4). Using the central value of n s and requiring w re <1/3 gives the bounds: r > 0.03 
(for p = 2 and p = 3) and r > 0.02 (p = 4). The region in parameter space that is associated 
with these more likely values of w re then allows for fairly small r values. 

8 Discusion and Conclusions 

Inflation includes a wide variety of models that give similar predictions for the fairly small 
number of available inflationary observables. The physics of reheating can provide an ad- 
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Figure 12: Plots of N re and T re , for hilltop inflation with p = 4 and for three different 
values of p. Shading is as in Fig. 2 


ditional opportunity to break this degeneracy. While CMB fluctuations themselves do not 
supply direct probes of the physics during the reheating era, the details of reheating affect 
the predictions for inflation (and vice versa) because they determine the nature of the cosmic 
thermal history after inflation (see e.g., [75-82]). Although we do not know exactly what 
occurred during reheating, we can make reasonable assumptions such as that the average 
equation of state during reheating was very likely between 0 and 1/3. This leads to indepen¬ 
dent constraints on observables like n s and r, that can then be tested against CMB data. 

One can parametrize our ignorance about reheating in terms of an equation of state, w re , 
a length N re (measured in terms of number of e-folds elapsed from the end of inflation), and 
a final temperature, T re . For any given inflationary model, one can write relations between 
the specific model parameters, the amplitude of the scalar power spectrum A s . the spectral 
index n s , and the reheating parameters (w re , N re and T re ). These relations are derived 
by accounting for the total expansion history between the time the observable CMB modes 
crossed outside the Hubble radius during inflation and the time of observation, and employ¬ 
ing a continuity equation for the energy density during the different cosmological epochs. 
We also assume that w re is constant. For single-field models the derivation is particularly 
straightforward. The main results are summarized in Eq. (2.11) for T re and w re (and in 
Eq. (2.10) for N re and w re ) as a function of inflationary parameters and observables. 
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n s 

Figure 13: Plot of the parameter space in the r vs. n s plane for the three hilltop models 
with p = 2,3, and 4. Shading is as in Fig. 3 



°0.93 0.94 0.95 0.96 0.97 0.98 0.99 1.00 


n s 

Figure 14: Plot of the parameter space in the vs. n s plane for the three hilltop models 
with p = 2,3, and 4. Shading is as in Fig. 3 


- 24 





























0.25 


0.20 



Figure 15: We recreate a version of the Planck Figure 12 in their recent Constraints on 
Inflation paper, using their 1 and 2 a TT, TE, EE + lowP constraints on n s and r [13], 
but plotting the parameter space for models such that there exists a reheating solution for 
0 < w re < Jj, as opposed to Planck’s choice of parameter space for which there is a solution 
with Nk between 50 and 60. Following the conventions in Planck’s version of the plot, the 
green line is (f > 3 , the black is </ 2 , the pink is 0 4 / 3 , the yellow c/), the red 0 2 / 3 , the orange 
Starobinsky/ Higgs model, the puple region is natural inflation, and the green region is the 
quartic hilltop model. 


We consider a broad range for the equation of state parameter, —1/3 < w re < 1, and the 
corresponding limits on CMB observables for different inflationary models. We notice that a 
0 2 potential would favor relatively large values of r: a reheating model with w re < 1 implies 
r > 0.11; to allow for a reheating model with w re <1/3 which is very probable, requires 
r > 0.14. Since it appears BICEP2’s signal is dust instead of primordial gravitational waves, 
it is difficult to reconcile q i> 2 inflation with the data. We also consider Starobinsky/Higgs 
inflation, natural inflation and the hilltop models. For Starobinsky and Higgs inflation, re¬ 
quiring w re <1/3 corresponds to r > 0.004. Because natural and hilltop inflation models 
have two free parameters, there are ranges of parameter space that can fit well the data for 
any value of the reheating parameter w re . For natural inflation, we find that Planck’s 2 a 
bound on n s favors a tensor-to-scalar ratio r > 0.05 for w re <1/3 (Fig. 9). For the same 
range of w re , the hilltop model, on the other hand, allows for smaller r values, specifically 
r > 0.02 for p = 2, r > 0.007 for p = 3, or r > 0.003 for p = 4 (Fig. 13). 

We show this parameter space in Fig. (15), where we recreate a version of the Planck 
Figure 12 in their recent Constraints on Inflation paper, using their 1 and 2 a TT, TE, EE 
+ lowP constraints on n s and r [13]. To get their model parameter space they impose 
between 50 and 60. Instead here, we don’t specify A/, but plot the parameter space for which 
there exists a reheating solution with 0 < w re < |. Constraining models in this way, using 
w re , is a nice model dependent but straightforward and well motivated way of representing 
the parameter space. 
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To conclude we find that considering broad, well-motivated physical constraints on the 
reheating equation of state indeed allows one to narrow the viable parameter space for in¬ 
flation models, offering an improvement over merely specifying whether or not an inflation 
model can reproduce the correct predictions at the pivot scale. These methods will become 
increasingly effective with future more precise CMB data. 

Note added: 

Our analysis was initially performed considering the Planck 2013 results for the scalar power 
spectrum parameters. Just after completion, but before submission, Planck released their 
2015 data, so we have updated our analysis using the new observational bounds on n s and 
A s . All presented results are now based on the Planck 2015 data. While completing the first 
version of this work, it was brought to our attention that a similar approach was carried 
out by [98]. Some of the results on natural inflation were reproduced in [98] and, where 
there is overlap, we find agreement if we consider the Planck 2013 bounds on the scalar 
power spectrum parameters. Furthermore, when this work was near completion, two papers 
concerning Higgs inflation and reheating were released [99, 100]. We find agreement with the 
results of [99] if we consider their pivot scale (k p = 0.002Mpc -1 as opposed to 0.05 Mpc~ x ). 
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